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bX)' We study general models for holographic superconductors in Einstein-Gauss-Bonnet gravity. We 

^ ; 

find that different values of Gauss-Bonnet correction term and model parameters can determine the 
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order of phase transitions and critical exponents of second-order phase transitions. Moreover we 
find that the size and strength of the conductivity coherence peak can be controlled. The ratios 
ujg/Tc for various model parameters have also been examined. 
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I. INTRODUCTION 



flO 



The anti-de Sitter/conformal field theories (AdS/CFT) correspondence [1|, |2| has provided a theoretical 
framework to describe the strongly coupled confornial field theories through a weakly coupled dual gravi- 
tational description. Recently stimulated from this correspondence, a remarkable connection between the 
condensed matter and the gravitational physics has been discovered, for reviews see 3|-|5|. It was suggested 
that the spontaneous U{1) symmetry breaking by bulk black holes can be used to construct gravitational duals 
of the transition from normal state to superconducting state in the boundary theory [(| . One can look at a 
(2+l)-dinicnsional superconductor and see its striking image involving a charged black hole with non-trivial 

"hair" in (3+l)-dimensions [7]. This investigation was further carried out beyond the probe limit by consid- 

I 
ering the back- reaction of the field on the spacetime [&\. Gravity models with the property of holographic 

superconductor have attracted considerable interest for their potential applications to the condensed matter 



physics, see for example 9l-|27l|. At the moment when the condensation occurs in the boundary CFT and in 
the gravitational counterpart a non-trivial hair for the black hole is triggered, there appears a phase transition. 
It was argued that this phase transition belongs to the second order [7]. The phenomenological signature of 



this phase transition was recently disclosed in the perturbation around such AdS black holes 



28|, 



Franco et al. recently introduced a generalization of the basic holographic superconductor model in which 



the spontaneous breaking of a global U(l) symmetry occurs via the Stiickelberg mechanism 30j. They found 
that a generalized Stiickelberg mechanism of symmetry breaking allows for a description of a wider class of 
phase transitions. This framework allows tuning the order of the phase transition which can accommodate 
the first order phase transition to occur, and for the second order phase transition it allows tuning the values 
of critical exponents 3l|. An interesting extension was done in 1321] by constructing general models for 
holographic superconductivity. It was found that except some universal model independent features, some 
important aspects of the quantum critical behavior strongly depend on the choice of couplings, such as the 
order of the phase transition and critical exponents of second-order phase transitions. In addition to the 
numerical investigation, analytical understanding on the phase transition of holographic superconductor was 



also provided in 



33|. 



It is of great interest to generalize the investigation on the phase transitio n in the holographic superconductor 
to the Einstein-Gauss-Bonnet gravity. It was observed in recent works 
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34| that the Gauss-Bonnet 



coupling affects the condensation and the higher curvature correction makes condensation harder to form. 
Further the high curvature correction also causes the behavior of the claimed universal ratio uj /Tc ~ 8 
unstable. In this work we are going to examine the effect of the Gauss-Bonnet correction on the order of the 
phase transition and critical exponents. We will also study the conductivity in a certain range of parameters 
and disclose the influence given by the Gauss-Bonnet coupling. 

The plan of the work is the following. In Sec. II, we will study the phase transitions of a general class of 
the holographic superconductor models via the Stiickelberg mechanism in the Gauss-Bonnet AdS black hole 
background. In Sec. Ill we will explore the effects of the Gauss-Bonnet correction terms and other model 
parameters on the conductivity. We will conclude in the last section of our main results. 

II. GENERAL SUPERCONDUCTING MODELS IN GAUSS-BONNET ADS BLACK HOLE 

We will consider the background solution of a neutral black hole in d dimensional Einstein-Gauss-Bonnet 
gravity [35|, |36 1 



ds' 



~f{r)dt' 



fir) 



r dxjdx^ 



(1) 



with 



^<^) - -. 



/ 4a / AfL2\ 



(2) 



where Ad is a constant of integration relating to the black hole horizon by r+ — (AfL^)^/^'*^^^ a is the Gauss- 
Bonnet coupling constant and L is the AdS radius. The Gauss-Bonnet correction a has an upper bound called 
the Chern-Simons limit a = L- /A, and a lower bound determined by the causality 
(1) goes back to the Schwarzschild AdS black hole. 
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In the limit a ^- 0, 



Consider a U{1) gauge field and the scalar field coupled via a generalized Stiickelberg Lagrangian 

S = I d'^Xy/^ 



m 



If^.F^'' ~ ^d^^d'^i^ - \ni^^ - \mW^.P - A^){d^p - A^) 



(3) 



with gauge symmetry A^ — > A^ + d^K and p— >p + A. 3^ is a general function of ij] which has the following 
form 



^{^) = ^2 ^ c^V> + Ci^^ 



(4) 



with the model parameters c^, 7 and C4. When c^ and C4 are zero, it reduces to the model considered in 

mm- 



Using the gauge freedom to fix p = and taking the ansatz ip = ip, At — (j> where ip, (j) are both real 
functions of r only, we can obtain the equations of motion 



f_ d-2 

J + ^-J- ■ 2/2" '^' f 



V'" + ( ^ + ^ ) ^' + ^S'W - — V' = , 



_ , 0. (5) 

r f 

These two equations can be solved numerically by doing integration from the horizon out to the infinity. At 

the asymptotic AdS boundary (r — )■ oo), the solutions behave like 

with 



A± = -[(d- 1) ± ^{d-lY+Am^Ll^ ] , (7) 



where L^g = 2a/ (1 — -^/l — Aa/ L^ ) is the effective asymptotic AdS scale 
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20[, ^ and p are interpreted 



as the chemical potential and charge density in the dual field theory respectively. Notice that both of the 
falloffs are normalizablc for f/;, so one can impose boundary condition that either ?/;+ or ^_ vanishes [7|,|8|. For 
simplicity, we will take ■(/;_ — 0. Moreover, we will set d—'a and m?L^ = —3 for concreteness. As a matter of 
fact, the other choices will not qualitatively modify our results. Thus, the scalar condensate is now described 
by the operator (0+) = ijj+ and we will discuss the condensate (0+) for fixed charge density. 

We will investigate how the phase transition depends on the coefficients c-y and C4. Due to the special 



interest in the case ^{ijj) = ijj^ + 04-0 '^ 
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18| , in (Ul we are going to set c-y = for the moment and pay more 
attention on the influence of C4 on the phase transition. Solving the equations of motion numerically, in Fig[T] 
we plot the condensate around the critical region for chosen values of C4 and different Gauss-Bonnet constants. 
For < C4 < 0.3, the transition is second order and the condensate approaches zero as (0+) ^ [Tc — T)^ , 
with mean field critical exponent (3 = 1/2 for all values of a. For C4 > 1, we observe that (0+) becomes 
multivalued near the critical temperature and the condensate does not drop to zero continuously at the 
critical temperature. This behavior keeps for all values of a. The analogous phenomenon holds as well when 
we consider the {0-) condensate but with different C4 range. In 31], for the case without Gauss-Bonnet 
constant, it was argued that the behavior for C4 > 1 indicates that the phase transition changes from the 
second order to the first order at C4 = 1. Here we find that the Gauss-Bonnet constant does not alter the 
result when C4 > 1. Choosing C4 £ [0.3, 1.0], we observe in Fig. [T]that for fixed C4, the transition point of the 
phase transition from the second order to the first order appears easily for the bigger value of a. In table U 



we list the critical value of ac separating the second order and the first order phase transitions for selected 
C4 within the range [0.3,1.0]. Since we concentrate on the five-dimensional spacetime, the separation point 
between the second and the first order phase transitions does not coincide with that in four dimensions, for 
example (04 = 0.9, ac = 0) here and (c4 = 1, a = 0) in '31i|. With the increase of C4, we see that ac becomes 
smaller. Thus we find that when C4 G [0.3, 1.0], not only C4 but also the Gauss-Bonnet constant can tune the 
order of the phase transition. The Gauss-Bonnet constant provides richer physics in the phase transition. 



TABLE I: The critical value of Qc which can separate the first- and second-order behavior for different 5('0) — ip'^ + Ciip'^ 
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FIG. 1: (Color online) The condensate < 0+ > as a function of temperature with fixed values C4 {c^ = 0) for different 
values of a, which shows that a different value of a can separate the first- and second-order behavior. The five lines 
in each panel from left to right correspond to decreasing a, i.e., 0.25 (black), 0.2 (red), 0.1 (green), (blue) and —0.1 
(orange). 



We also have the interest to see the influence of c^ on the phase transition. We will concentrate on 7 in 
the range 3 < 7 < 4. At this moment we set C4 = 0. In Fig. [5] we exhibit the condensate of < 0+ > for 
selected values of Cy, a and changing 7. We see that for fixed small values of a and Cy, it is more possible to 
observe the appearance of the first order phase transition when 7 is smaller. The a influence on the order of 
the phase transition is also exhibited. Bigger a will cause the first order phase transition to appear easier for 



fixed c^ and 7. This is consistent with the case for nonzero C4 but zero c^. The effect of c^ is similar to the 



Gauss-Bonnet constant, bigger c-^ brings the first order phase transition easier for fixed 7 and a. 
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FIG. 2: (Color online) The condensate < 0+ > as a function of temperature with fixed values c^ and a for different 
values of 7, which shows that the type of phase transition depends on the Gauss-Bonnet correction term a, model 
parameters 7 and c~f. The three lines in each panel from left to right correspond to decreasing 7, i.e., 4.0 (black), 3.5 
(red) and 3.0 (blue and dashed). 



In addition to showing that the gravity duals can lead to both first and second order phase transitions, it is 
also of great interest to examine in the second order phase transition whether different choices of 0-^,04,7, a 
result in different critical exponents from the prediction of the mean field. When c^ — ^ and C4 = 0, the 
critical exponent was shown in agreement with the mean field value /3 = 1/2 in the Einstein-Gauss-Bonnet 



gravity 



m 



20|. An interesting behavior arises for c-y < 0. In Fig. [31 for 



-1 and C4 — 1/2, we present 



the condensate < 0+ > as a function of 1 — T jTc in logarithmic scale with different values of a for choosing 
7 = 3.0, 3.25, 3.5 and 4.0 respectively. Three lines in each panel from the bottom to the top correspond to 
a = 0.25, —0.1 and 0.1. We see that the slope is independent of the Gauss-Bonnet correction term a but 
sensitive to the model parameter 7. Further analysis shows that, near the critical temperature Tc, the critical 
exponent /3 ~ 1.00 for 7 = 3.0, l3 ~ 0.79 for 7 = 3.25, (3 ~ 0.67 for 7 = 3.5 and l3 ~ 0.50 for 7 = 4.0, which is 



independent of a! The relation between the critical exponent /? and the parameter 7 can be expressed as 

1 



P 



7-2 
This behavior is consistent with that seen for the AdS4 black hole with the scalar mass m^L^ 



31 



(8) 



which tells us that the critical exponent /3 depends only on the model parameter 7 but is independent of the 
scalar mass and the background spacetimes. 
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FIG. 3: (color online) The condensate < 0+ > vs 1 — T/Tc in logarithmic scale with different values of a for 7 = 3.0 
(top and left), 3.25 (top and right), 3.5 (bottom and left) and 4.0 (bottom and right). The three lines in each panel 
from bottom to top correspond to a = 0.25 (blue), —0.1 (red) and 0.1 (green). These panels show that the slope is 
independent of a but sensitive to 7. 



In Fig. |4]we present the condensate as a function of temperature for c.^ = — 1 and C4 = 1/2. It exhibits 
that the critical temperature Tc is independent of the model parameter 7 but depends on the Gauss-Bonnet 
correction term a, i.e., T^ = 0.209 for a = -0.1, Tc = 0.198 for a = 0, Tc = 0.186 for a = 0.1 and Tc = 0.158 
for a = 0.25, which shows that the positive Gauss-Bonnet correction will suppress the condensation but the 
negative one will enhance it 



2l|. In fact, we note that the critical temperature Tc is not sensitive to the 
coefficients c^ and C4 in the general function di'ip) from Figs. [1] [5] and |3] for the phase transition of the second 
order. Thus, we point out that the critical temperature Tc depends on the background spacetimes but not on 



the model of g', 
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FIG. 4: (color online) The condensate < 0+ > as a function of temperature with fixed values a for different model 
parameters 7. The four lines in each panel from bottom to top correspond to increasing 7, i.e., 3.0 (blue and dashed), 
3.25 (green), 3.5 (red) and 4.0 (black). These panels show that the critical temperature Tc is independent of 7 but 
depends on a. 



III. CONDUCTIVITY 



Now we jump to investigate the influence of the Gauss-Bonnet correction term a, model parameters 7, c-y 
and C4 in ^{iJj) on the conductivity. 

In order to calculate the conductivity, we consider the perturbed Maxwell field SA^ = Ax{r)e''^'^*dx. The 
equation of motion for SA^ reads 



a: 



f r 



AL 



P f 



A^=0. 



(9) 



We still restrict our study to d = 5 in order to avoid the complicated behavior in the gauge field falloff 
in dimensions higher than five. We solve the above equation by imposing the ingoing boundary condition 
Axir) ^ firf^ near the horizon. In the asymptotic AdS region, A.^ = yl(o) + 4^ + d!!!f^lLk 121^1 j^ 
should be noted that the appearance of the arbitrary integration constant. A, leads to a logarithmic divergence 
in the retarded Green's function G^ which gives the conductivity a. In order to remove the divergent term 



\m, 



logAr, we add a boundary counter term in the gravity action as suggested in _395, which can specify the 
renormahzation scale when regulating the action [9|. With this 



appropriate boundary counter term to cancel 



the logarithmic divergence, we can express the conductivity as [9, 
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20| 



lU! 



(10) 



For the general forms of function ^{i/j) = i/j'^ + c-^tp'' (setting C4 = for clarity), one can obtain the conductivity 
by solving the Maxwell equation numerically. We will focus on the case for the fixed scalar mass m?L^ = —3 
in our discussion. 
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FIG. 5: (color online) The real part of the conductivity for (3 + l)-diniensionaI Gauss-Bonnet superconductors with 
different values of 7, c^ and a. 



In the top-three panels in Fig. [5] we have plotted Rea{ijj) by fixing 7 = 4 for different values of c^ and 
a. It clearly shows that for fixed a the coherence peak gradually becomes stronger and narrower with the 
increase of c^. This is consistent with the result of C4 observed in 31], which indicates that c^ controls the 
magnitude of the fluctuations of the condensate. For bigger c^ we have already seen that the fluctuations 
are strong enough to induce the first order phase transition. In the bottom-three panels in Fig. [5] we show 
the result by fixing c^ = 1 and varying 7 with selected a at T jT,^ ~ 0.7. We observe that with the increase 
of 7 the coherence peak increases for the same a, which shows that the model parameter 7 also controls the 
magnitude of the fluctuations. From Fig. [5l we also learn that for the selected c-^ and 7, the coherence peak 
becomes lower with the larger a, which shows that the fluctuations in holographic superconductors will be 
suppressed by the higher order curvature corrections. 

In 9| it was argued that there is a universal relation between the gap Wg in the frequency dependent 



10 



conductivity and the critical temperature T^'. ^g/T^ ~ 8 respected to a good approximation by all cases 
considered. However this claimed universal relation was challenged when the higher curvature corrections are 
19|, |20f. The available discussions were in the model when ^{ip) — ip^. In the general ^ 



taken into account 



form in our work, we will show that the model parameters 7, c^ also modify the claimed universal relation 
even when a — 0. In Fig. IHlwc plot the conductivity at temperature T/Tc ~ 0.3. The blue (bottom) line 
represents the real part of the conductivity and the red (top) line is the imaginary part. We can easily find a 
gap in the conductivity with the gap frequency ujg changes with the values of Gauss-Bonnet correction term 
a, model parameters 7 and c-y. Fixing 7 and c^, the gap frequency ujg becomes larger for bigger a, which 



agrees with the finding in 



19 



201 . The gap frequency cug increases with 7 for fixed c^ , a and grows with c^ for 



selected 7, a. The deviation from cOg/Tc = 8 becomes bigger with the increase of 7 and Cj. This holds even 
when a = 0. This shows that not only the high curvature correction, but also the form of the scalar field ^ 
will affect the so-called universal relation ujg/Tc ~ 8 [9|. 
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FIG. 6: (color online) Conductivity for (3 -I- l)-diniensional Gauss-Bonnet superconductors with fixed values of a. 
for different models with 5'(V') = V"^ + c^i^'^ ■ The blue (bottom) line and red (top) line represent the real part and 
imaginary part of the conductivity respectively. 



11 

IV. CONCLUSIONS 

We have introduced a general class of gravity dual with Gauss-Bonnet corrections to describe both first and 
second order phase transitions at finite temperature in strongly interacting systems. In the probe limit, we 
found that besides the parameters which define 5' can separate the first and second order phase transitions, 
different values of the Gauss-Bonnet correction bring richer physics and can also change the order of the phase 
transition in the generalized system. However for the second order phase transition, we observed that the 
shift of the critical exponents from that of the mean field result only appears for the parameters defining ^ 
and is independent of the Gauss-Bonnet coupling. This is different from what we observed for the critical 
temperature Tc, which only depends on the Gauss-Bonnet constant while has nothing to do with other model 
parameters. 

We also discussed the influences of the Gauss-Bonnet corrections and other model parameters on the con- 
ductivity. We found that the size and strength of the coherence peak can not only be controlled by the 
parameters that define ^ as observed in [SI], but also be infiuenced by the Gauss-Bonnet coupling. This 
shows that ^ together with the high curvature correction controls the magnitude of the fluctuations in the 
system. Furthermore we examined the relation between the gap cjg in the frequency dependent conductivity 
„a .e e„«c. _.„ ... we .„.. .H. . ..... . .He H,„ e™.„e .„,.„e „. .He ...o .,r. 

mm, 



observed in 

UJg/Tc « 8. 



the model parameters defining ^ also give corrections to the so-called universal relation 
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